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This paper examines the performance of a new Marnoch heat engine, which uses dry air and a pneumatic 
piston assembly to convert thermal energy to electricity. The system has unique capabilities of operating 
over temperature differentials less than 100 K. Unlike a common Stirling engine, the heat exchangers and 
piston assembly are not co-located, which is beneficial for positioning of heat exchangers in various 
configurations. This paper presents an operational laboratory-scale, proof-of-concept Marnoch heat 
engine (MHE), including its performance and power generation capabilities. It also presents a thermo¬ 
dynamic analysis of the system. Based on the MHE results, component modifications are made to 
improve its performance. The configuration has an efficiency of about thirty percent of a Carnot heat 
engine operating in the temperature range between 272 1< and 372 K. Experimental data is acquired to 
provide verification of the predictive model, as well as demonstration of the MHE’s capabilities for 
efficient generation of electricity from waste heat sources. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

Rapid population growth in the world, at about 2% per annum, 
will create a large growing need for energy [1,2]. World energy 
demand in the past decade has increased by 20% in comparison to 
the 1990s [2]. Total daily world oil demand is predicted to increase 
from 77 million barrels in 2003 to about 119 million barrels in 2025, 
a 55% increase [3], In 2004, a research study by the Pembina 
Institute in Canada [4] showed that electricity consumption during 
peak hours by the year 2020 could be reduced by 30% below 2004 
levels, through a series of energy efficiency and demand reduction 
policies. The International Energy Agency (IEA) reported that more 
than 60% of the total worldwide energy demand is being supplied 
by fossil fuels [5], In the U.S., the fossil fuel dependence was about 
85% in 2007 [6], Utilizing systems that work on more efficient 
thermodynamic cycles reduces the amount of pollution emitted to 
the environment [7]. Waste heat recovery from thermal power 
plants, industrial facilities, furnaces, and so forth, represents an 
important element of worldwide efforts to reduce greenhouse gas 
emissions and fossil fuel dependence [8]. Various types of heat 
engines are capable of recovering heat from lost energy. This paper 
focuses on heat recovery from a new type of heat engine, called 
a Marnoch heat engine (MHE). 
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Most gas power cycles were originally developed in the nine¬ 
teenth century. Carnot and Stirling cycles are two well known power 
cycles [9]. In a Carnot cycle, the system undergoes a sequence of four 
reversible processes. This includes two adiabatic processes and two 
isothermal processes. The Stirling cycle represents the processes in 
a heat engine with regenerators (see Fig. 1). In this cycle, compres¬ 
sions and expansions are isothermal. Heat absorption and rejection 
in a Stirling cycle occur through constant volume processes [10]. In 
comparison with ICEs (internal combustion engines), heat engines 
have the advantages of operating on any heat source, not solely 
combustible fuels, simpler mechanisms, among others. 

In Stirling engines, expansion and compression takes place in 
a cylinder with a power piston. A displacer transfers the working 
fluid through the heater, cooler, and regenerator. Fig. 2a shows the 
Alpha configuration of a Stirling engine. This assembly has no dis¬ 
placer, and the two pistons are working together so that volume 
stays constant. Cold and hot pistons move uniformly in the same 
direction and they provide constant volume heating/cooling of the 
working fluid. When the majority of gas is in one cylinder, the other 
piston expands and compresses the working fluid. Expansion 
occurs on the hot side, and compression on the cold side [11,12]. In 
the Beta configuration (Fig. 2b), a displacer and power piston are 
installed within the same cylinder [13], The working fluid between 
the heat source, sink, and regenerator is moved by the displacer. 
The Gamma configuration (Fig. 2c) has two separated cylinders: 
one for the displacer and the other for the power piston [14]. In 
these configurations, the Stirling engine suffers from several 
disadvantages, including low specific power (large size with 
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Nomenclature 

Cs 

cold source 



d 

destroyed 

A 

area [m 2 ] 

e 

electric 

AD 

angular distance [o] 

f 

final 

B 

bias limit 

fa 

from air 

Cave 

constant 

fpc 

friction between sealing and the cylinder body 

Ex 

exergy [/] 

ftp 

friction between rack and pinion 

Ex 

exergy rate [W] 

fw 

from water 

h 

specific enthalpy [//mo/[ 

gen 

generated 

H 

enthalpy [/] 

Hs 

heat source 

Jm 

rotational mass moment of inertia [ kgm 2 ] 

i 

point number 

k 

radius of gyration [m] 

in 

inlet 

KE 

kinetic energy [/] 

int 

initial 

L 

length [m] 

loss 

loss 

Q 

heat flow rate [W] 

out 

outlet 

M 

mean value 

p 

piston 

m 

mass [kg] 

Q 

heat 

m 

mass flow rate [kg/s ] 

th 

thermal 

N 

number of experiments 

to 

total 

P 

precision limit 

ta 

to air 

R 

universal gas constant (J/moM<) 

tw 

to water 

S 

entropy [ l/kg K] 

w 

water 

T 

temperature [/C] 



t 

time (s) 

Greek letters 

V 

volume [m 3 ] 

r 

exergetic temperature factor 

u 

uncertainty 

T 

torque (N.m) 

w 

work output [/] 

V 

efficiency 

X 

piston length [m] 

A 

change 



a 

standard deviation 

Subscripts 



0 

dead state (ambient condition) 

Acronyms 

1, 2 ,... 

state 1, 2 , ... 

NSP 

number of strokes per minute 

n 

second law 

Rt 

rotational time [s] 

a 

air 

RPM 

revolutions per minute 

apl 

all piston losses 

SRDR 

shaft rotation angle per revolution [°] 

ave 

average 




relative small power output), requires large temperature differen¬ 
tials to operate efficiently, and so on. 

Various other types of low-temperature differential heat engines 
(LTDs) have been developed previously to generate mechanical 
power from low-temperature differentials. Abdullah et al. [15] 
developed a solar-based LTD, which generated power of about 
100 W [15]. Kongtragool et al. [16] reported the characteristics of 
different LTD engines. It has been shown that the torque output of 
LTDs has been generally insufficient to drive industrial electric 
generators with high torque [17—19]. This paper presents a new 
type of heat engine, called a Marnoch heat engine (MHE), which has 
promising potential to overcome some limitations of past LTD 
technologies. It is capable of producing electricity from low grade 
heat sources and temperature differentials, below 100 K. Operation 
over such small temperature differences allows for a highly flexible 
system. This paper will acquire new experimental data and develop 
predictive models of the MHE performance, in order to demonstrate 
its promising performance as a new type of LTD for waste heat 
recovery. 

2. Description of Marnoch heat engine 

The Marnoch heat engine (MHE) operates similarly to other heat 
engines, by converting thermal energy to mechanical work, fol¬ 
lowed by electricity generation. The MHE consists of four cylin¬ 
drical heat exchangers, each fitted with helical tube coil heat 


exchangers. Each heat exchanger is connected to both a hot source 
and a cold sink, via tubes and valves that allow a controlled mass 
flow at any step of the process. Initially, all shells are kept at 
a constant pressure and mass. A half of the Marnoch cycle is driven 
by heat supplied to the exchangers, which drive pistons external to 
the exchangers. The other half of the cycle connects to a colder 
exchanger, facilitating thermal expansion of the gas in the hot 
exchanger and thermal contraction in the cold exchanger, via 
a piston assembly. The heated gas flows through the piston 
assembly of the heat engine and produces shaft work. 

The configuration of the MHE with two pairs of heat exchangers 
is shown in Fig. 3. This figure shows the current assembly of the 
MHE. Fig. 4 indicates how the heat exchangers are connected 
internally and externally to the heat source and sink. When oper¬ 
ation starts, hot water moves from the heat source through heat 
exchanger 1 and heats the gas (air) inside the shell. The air pressure 
inside the shell increases due to the temperature increase in the 
vessel. The coolant flows from the heat sink through to tubes of 
heat exchanger 2. As the result, it reduces the air pressure in the 
heat sink. 

In the next step, when the air temperature inside the shells 
reaches a desired level, the valves are opened. Consequently, the air 
flows from the high pressure to low-pressure shell. Pipes are con¬ 
nected for a piston-cylinder assembly in series. Therefore, when air 
flows from the heated heat exchanger to the cooled heat exchanger, 
it moves the piston. When the piston reaches its end position, the 
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Fig. 3. Experimental apparatus of the Marnoch heat engine. 


Fig. 1. Stirling cycle (a) P-V diagram; (b) T-s diagram. 
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Fig. 2. Different configurations of Stirling engines. 


pneumatic valves redirect high pressure air to the other side of the 
piston to move it back to its initial position (see Fig. 4). The air 
moves from the high pressure shell to the heat sink, and it is stored 
in the cold shell. 

By transferring air from the high pressure side and storing it in 
the cold heat exchanger, the pressure of the cold side increases 





























































































































6882 


P. Saneipoor et al. / Energy 36 (2011) 6879-6889 


gradually, and a larger mass of air accumulates on the cold side. 
After a certain number of strokes, based on the size of the 
exchangers and pre-charged pressure, the pressures of cold and hot 
sides become nearly equal. The air pressure is no longer sufficient to 
move the piston, after which a PLC (programmable logic controller) 
disconnects the first pair of exchangers and connects exchangers 3 
and 4, for the remainder of the cycle. 

In the final operation step, the roles of heat exchangers 1 and 2 
change. During this period, the engine is using heat exchangers 3 
and 4 to generate the pressure difference, for the engine to operate 
continuously. Throughout the second mechanical cycle of heat 
exchangers 1 and 2, heat exchanger 1 is used as the heat sink, and 
heat exchanger 2 is used as the heat source. Solenoid valves are 
changing the direction of water from the heat source and sink, such 
that hot water flows to heat exchanger 2 and cold water flows to 
heat exchanger 1 during the cycle. When the pressure inside heat 
exchangers 3 and 4 is nearly equalized, the first pair is reused to 
generate the temperature difference, and the last step is repeated 
for the second pair. 

Thermodynamic and mechanical considerations show that there 
are major differences between the Stirling engine and MHE. In Stir¬ 
ling engines, a part of the mechanical work output from the flywheel 
is used in the compression process after each stroke. However, in the 
MHE, all of the energy output of the piston expansion is transferred to 
the flywheel. The reduced energy output due to the back pressure on 
the piston is equal to the compression energy to complete the ther¬ 
modynamic cycle. Also, in MHEs, piston-cylinder assembly and heat 
exchangers can be located away from each other, whereas in the 
Stirling engine, the heat exchangers and piston assembly are 
attached. This feature of the MHE results in its ability to use a variety 
of sources to generate pressure differences from temperature 
differentials. Thirdly, unlike Stirling engines, after each stroke in 
MHEs, the enclosed area of a P—V diagram changes. This area is 
reduced, because after each stroke, the expansion pressure reduces 
and the compression pressure increases. Fourthly, in the mechanical 
configuration of the MHE, no regenerator is required. These differ¬ 
ences will lead to noteworthy advantages of the MHE over a Stirling 
engine. 

Fig. 7 demonstrates the sequential processes of the MHE cycle. 
Process 1—2: dry air inside the shell absorbs heat from the liquid 
flowing through the tube bundles. The shell is a rigid body, thus, the 
air volume stays constant and the gas pressure increases. Process 

2— 3: in this process, the valve that connects the high pressure shell 
to the low-pressure shell, via the piston assembly, is opened. Air 
expands and the piston moves forward. The expansion is almost 
isothermal, because the hot liquid flows through the tube bundle to 
keep the gas temperature constant during the expansion. Process 

3— 4: dry air accumulates in the low-pressure shell and it is cooled 
with the tube bundles that transfer cold water to the heat 
exchanger. The volume of the gas inside the shell stays constant and 
the pressure decreases due to the gas temperature change. Process 

4— 1: dry air is compressed with the piston and simultaneously 
cooled down with the liquid that flows through the tube bundles. 
This compression process is isothermal. 

3. Experimental apparatus of the Marnoch heat engine 

A proof-of-principle prototype was constructed, using a pinion 
rotary actuator assembly. The body of the cylinder is made of 
aluminum, which is hard-coat anodized and permanently sealed 
for wear resistance and endurance. The sealing in the system is 
critical to avoid leakages from the piston assembly to the 
surroundings, and from one side of the piston to the other side. 
Also, the friction loss between the piston and cylinder must be 
minimized. In the current experimental prototype, the maximum 



Direction of piston motion 

Fig. 5. Schematic of the MHE with one pair of exchangers for the thermodynamic 
analyses; section (a) is for the charging period and section (b) is for the discharging 
(operation) period. 


operating pressure is 1723.7 kPa, with a breakaway pressure of 
34.5 kPa. 

Rotation of the shaft by each complete stroke is 270°. The 
displacement per degree of rotation in the cylinder is 1.23 cm 3 . The 
piston in this prototype has a diameter of 63.30 mm. Therefore, in 


3 Carnot Cycle 



Marnoch Cycle 
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Fig. 6. Comparison of Carnot and Marnoch cycles. 
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Operation Cycle 1 Operation Cycle 2 


Isovolumetric 



Operation Cycle 3 Operation Cycle 4 

Fig. 7. Four power cycles of the MHE for five sequential piston revolutions. 


each stroke, 331.83 cm 3 of the working fluid is transferred to the 
low-pressure shell. Seals in the piston are designed to continually 
operate between temperatures ranging from 255 I< to 355 K [20], 

The torque outputs of the piston-cylinder, in four pressure 
differences, are shown in Table 1 for two different diameters. These 
values refer to a piston-cylinder with a diameter of 63.30 mm, in 
the current MHE, and they will be compared with the torque output 
from a piston assembly with a diameter of 81.21 mm. In Table 1, the 
higher torque output of the system requires higher pressure 
differences. For piston sizing of the MHE, if the volume of the 
cylinder is large, consequently the amount of working fluid trans¬ 
ferred from one exchanger to the other is high, per stroke. Hence, if 
the volume of the heat exchangers does not match the piston size 
after a few strokes, the engine operation will stop and another pair 
of exchangers must be charged for operational use. The charging 
mode is time consuming, and delays in piston operation happen 
when exchangers are switched. When a larger piston-cylinder is 
used, less time is required for changing the direction of the piston 
through the valves. Thus, the power output from the piston 
assembly will be higher. Moreover, the torque output is more than 
the output of smaller pistons. When smaller pistons are used, the 
direction must be changed more frequently. Changing the direction 
of valves is time consuming, so a specified amount of electricity is 
produced over a longer duration. Therefore, the power output of 
the piston is lower. 

The inaccuracies related to measurement devices are bias errors. 
The measurement accuracy of the pressure gauges and power 
meter is about ±1%. Random variations are errors that occur due to 
uncontrolled parameters when making repeated measurements. 
They are related to the precision of the experimental readings and 


Table 1 

Torque output of the piston [20], 


Pressure difference (kPa) 

344.74 

517.11 

689.475 

1723.69 

Torque output (N.m) 

(Piston diameter: 63.30 mm) 

24.29 

36.38 

48.58 

121.35 

Torque output (N.m) 

(Piston diameter: 81.21 mm) 

64.40 

96.04 

128.92 

322.23 


determined by finding the deviation of readings. The root-sum- 
square method was applied to find the measurement uncer¬ 
tainties. Using the method of Kline and McClintock [21], the 
measurement uncertainty (U) for the experimental data is calcu¬ 
lated by the following equation: 

U = Vb 2 + P 2 (1) 

where B is the bias limit, and P is the precision limit for the 
experimental readings. Bias is the error associated with the inac¬ 
curacy of the measuring equipment. The precision can be found 
based on the following equation [22]: 


where a and N are the standard deviation and the number of 
experiments, respectively. The measurement uncertainties will be 
reported in the next section, alongside experimental results that 
are presented in the figures. The estimated error in the experi¬ 
mental data is approximately ±4.6%. This will place the experiment 
within a 95% confidence interval (see Table 2). 

4. Thermodynamic analysis 

Both energy and exergy analyses are performed to evaluate the 
performance of the MHE [23]. The schematic of the system is shown 
in Fig. 5, which represents the first operation cycle in two different 


Table 2 

Estimated uncertainty of the measured power output from the generator and 
actuator. 



Precision 
error (W) 

Bias 

error (W) 

Uncertainty 

(W) 

Uncertainty 

(%) 

Total power output 
from generator 

48.5 

19 

52.1 

4.3 

Total power output 
from actuator 

114.9 

47 

124.2 

4.0 

Peak power output 
from generator 

3.6 

1.0 

3.8 

4.6 
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stages (charging and discharging). In Fig. 5a, both air valves are 
closed and initially both shells have the same pressure. When water 
flows through the heat exchangers, heat exchanger 1 is heated to 
raise the pressure, and at the same time, heat exchanger 2 is cooled 
and its pressure decreases. In the second step (Fig. 5b), the air 
temperatures inside both heat exchangers are close to the temper¬ 
atures of the water streams. Air valves are opened and the set of heat 
exchangers and piston start their operation. Operation of this set 
continues until the air pressures inside both heat exchangers 
equalize. During charging and discharging periods, water flows 
through the heat exchangers continuously, and it does not shut off 
until the last stroke. This occurs to achieve near isothermal expan¬ 
sion and compression strokes of the piston. 

4.3. Mass and energy analyses 

Heat exchanger 1 is a shell and tube heat exchanger with both 
fluids unmixed. A heating fluid (hot water) flows through the tubes 
and exchanges heat with the working fluid (air) in the shell side. 
The overall mass balance can be written as 

Am = ™in ~ ™ out ( 3 ) 

Eq. (3) is written separately for water as 

m w i = m w2 (4) 

For air, the system is not at steady state, and the conservation of 
mass requires 

^+^5 = 0 (5) 

From Eq. (5) the total amount of air mass that leaves the hot heat 
exchanger can be calculated (see Fig. 5b). Integrating over time 
gives the change of mass in heat exchanger 1 during the overall 
process: 

t 

J( m ^y t = m initial -m flnal (6) 

o 

In Eq. (4), mj mna / denotes the initial air mass in shell 1, and mj\ na i 
denotes the air mass that remains in shell 1 after one cycle of 
operation. For the whole operation, over a period of time (t), the 
continuity equation for the transient process is written as 

(m initial ~ m fina i ) = - XI m 5 (At) (7) 

In the first cycle of operation, as shown in Fig. 5, heat exchangers 
1 and 2 are connected to the hot and cold reservoir, respectively. 
The general energy balance is given by: 

AE = tin-J2 tout (8) 

where AE is the total energy change in the shell and tube heat 
exchanger. This includes the energy change in both liquid and gas 
states. Also, J] E in is the summation of the total energy input to the 
heat exchanger, and J2^out is the energy that leaves the heat 
exchanger control volume. This equation is used for heat 
exchangers 1 and 2 for both modes: charging and discharging. In 
the first step for heat exchanger 1, heat is transferred into the 
system by the hot water stream. A portion of its energy is trans¬ 
ferred to the gas inside the heat exchanger. A fraction of J2 £)'n g° es 
back to the surroundings by the exit water stream, and the 
remainder is lost to the surroundings. Thus, the energy balance for 
the water stream becomes 


lii w ] h w ] — m w2 h w2 + Q at + Q/oss (9) 

where the values of rii wl and rh w2 are obtained from Eq. (4), and h w i 
and h w2 denote the enthalpy of water that enters and exits the heat 
exchanger, respectively. The term m w2 h w2 (energy flow to the 
ambient via the outlet water stream) changes over time and its 
value depends upon Q_ ta (heat transfer rate from the water to the 
working fluid). This rate decreases as the temperature of the gas 
approaches the water temperature. Also Q (oss is the amount of heat 
rejected to the ambient from the hot stream tubes, and it also varies 
with time because the fluid temperature changes with time. The 
heat transfer rate (Q ta ) from the water stream to the air is 
expressed as 

t 

J Q ta (dt) = mc v AT (10) 

o 

When Eq. (10) is used for the charging period, the air mass (m) is 
constant inside the shell, and its value is in Eq. (7). 

When the piston operation starts, in Fig. 6, it is assumed that 
Ti = T 2 = T 5 , where Tj and T 2 are the water inlet and outlet 
temperatures respectively, and T5 is the gas temperature inside the 
hot heat exchanger. The heat loss from the pipe, which connects 
shell 1 to the piston assembly, is neglected (well insulated system). 
Hence, air temperatures at point 5 (Fig. 5) and inside shell 1 are the 
same. 

When air leaves the hot heat exchanger, the pressure of the shell 
decreases. The expansion process in this heat engine is kept as close 
as possible to isothermal to obtain maximum work output. Water 
flows through the heat exchangers during piston operation, so 
when the gas is expanded, it is being heated simultaneously. The 
energy balance equation inside exchanger 1 during the operation 
mode can be expressed as 

AE = Q ta -m 5 h 5 (11) 

where AE is the rate of energy change in exchanger 1 due to heat 
and mass transfer. 

4.2. Exergy analysis 

To determine the performance of the system at different 
ambient conditions, an exergy analysis is also performed. The 
physical exergy of each point Ex,- can be calculated by 

Ex, = (H f - H 0 ) - T 0 (S ; - S 0 ) (12) 

where H 0 , T 0 and So are enthalpy, temperature and entropy of air 
under ambient conditions, respectively. Also H, and S, are the 
enthalpy and entropy at each point within the system boundary. 
The thermal exergy content is given by: 

Ex a = j (l- T T °)dQ (13) 

int 

where dQ is the incremental heat transfer, and the integral is taken 
from the initial state (inf) to the final state (J). This thermal exergy is 
the minimum work required in bringing the system to its final state 
from the dead state. In the MHE prototype, under most conditions, 
heating the exchanger starts from the dead state. 

T 

The dimensionless quantity (1-^) is called the “exergetic 
temperature factor”, and it is denoted as [24]: 

r=( l-y) (14) 
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For heat transfer across a segment (r) of the control surface, for 
which shell 1 is charged initially, the thermal exergy associated 
with heat transfer is 

Exq = J qr-rdAr (15) 

r 

where q r is the heat flow per unit area at a region on the control 
surface, and the temperature is denoted as T r . The dimensions of 
the copper tubes inside the heat exchanger are used to find the 
surface area (see Table 3). If the temperature (T) of the control 
mass is constant, the thermal exergy transfer associated with heat 
transfer is 

Ex a = (l = IQta (16) 

where To is the ambient temperature and T is the working fluid 
temperature during the charging mode. 

4.3. Shaft work output of the piston assembly 


v f vf 

W = J ^dV = C J IdV (19) 

Vint Vint 



The piston surface is constant during operation, but the length 
changes. Using AV = A AL in Eq. (20), it results in 


W = Cln(\ + ^J = C/n(l +C a „ e AL) (21) 

where C ave is a constant that is defined as C ave = y.. Each piston 
assembly in the MHE has a specific value of C ave . 

Eq. (20) shows the work output of the piston after one cycle. 
Then some air mass is transferred to the other exchanger. When the 
work output of the next stroke of the piston is calculated by Eq. 
(20), the work performed by the compression of the piston is 
negative. 


In Fig. 6, the thermodynamic power cycle of the MHE is shown 
and compared with an ideal Carnot cycle. Unlike the Stirling cycle in 
Fig. 6b, the enclosed area of the MHE is reduced after each stroke 
during the operation of one exchanger pair. Fig. 7 is a demonstra¬ 
tion of the area reduction of cyclic work output from the piston 
assembly. Pre-charge pressures and available temperature differ¬ 
ences will affect the number of piston cycles, and number of power 
cycles per operation period. The compression and expansion 
processes in the MHE are assumed isothermal, and air is assumed to 
be an ideal gas. The ideal gas law (PV = nRT) will be assumed to find 
the air properties. In this equation, P is the absolute pressure of the 
air, n is number of moles of air and R is the universal gas constant 
which is 8.314 (J/mol •K ). 

The MHE has a closed cycle, so no gas enters or leaves the system 
during the operation cycle. Hence, the number of moles of air (n) is 
constant. Also, the expansions and compressions are nearly 
isothermal. Therefore, the value of nRT is nearly constant, which is 
shown by C in Eq. (17). Thus, the ideal gas equation implies: 


P = 


C 

V 


(17) 


The work output of the piston as the result of one expansion can 
be found as 


Vf 

W = f PdV (18) 

Vint 

where P is a function of the volume. 

Substituting Eq. (17) into (18) yields 


4.4. Analysis of transmission system 

Table 2 shows the torque output from the current piston 
assembly at certain pressure differentials, based on the manufac¬ 
turer’s specifications [20]. To find the rotational velocity of the 
piston, the following equation is used: 

RPM P = (22) 

where RPM p is the piston shaft rotation per minute, NSP is the 
number of strokes per minute and SRDR is the shaft rotation 
degrees per revolution. In the current piston assembly, the pinion 
rotates 270° after each revolution. 

To find the power output for the piston, RPMp is multiplied by 
the torque output from the piston assembly, 

Power p = xp-(RPM p ) (23) 

where t p is the torque output of the piston assembly. The value t p 
changes after each revolution, because the pressure differences 
after each stroke decrease. In order to estimate the power output 
from the piston assembly, the output kinetic energy from the piston 
to the flywheel is calculated. Determining the losses in chains and 
bearings, the energy output from the piston can then be estimated. 

The thickness of the disc is negligible in comparison to its 
diameter. Thus, the kinetic energy from the flywheel is written as 

KE = \ j m w 2 (24) 

where I<E is the kinetic energy of the flywheel, J m is the rotational 
mass moment of inertia, and oj is the angular velocity. As given 
below: 


Table 3 

Helical coil specifications. 

Thickness of tube 

Inner diameter before deformation 

Outer diameter before deformation 

Number of coils 

Length of heat exchanger coils 

Radius of bends 

Total length of tubes 

Total occupied volume by tubes 

Percentage of occupied shell volume by tubes 


0.64 mm 
10.92 mm 
12.2 mm 
37 

630 mm 
65 mm 
17,350 mm 
2006.23 x 10- 6 mm 
2.57% 


W li 2 
~ 8 2 " 


(25) 


where w is the weight load, g is the gravitational constant, and k is 
the radius of gyration. The value of to to be substituted in Eq. (24) is 
calculated by 


AD 

to = — ■ (0.035) 


(26) 


where AD is the angular distance and Rt shows the rotational time. 
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4.5. Efficiencies of the Marnoch heat engine 


The thermal efficiency of the system is given by 


_ Woutput 
Qhs 


(27) 


where W ou tput is the work output of the cycle, and jjthis the thermal 
efficiency of the system. The variable W output can be calculated 
using Eq. (19). The value Qhs is the heat flow into the system from 
the heat source. It is defined as the total heat flow into the system 
from the hot water stream, which can be calculated by substituting 
Qta (amount of heat transferred to the air) into Eq. (27). Also, the 
overall efficiency of the system ( r \ t0 ) can be expressed as 


W e 

Qh 


where W e is the electric work output from the engine. 
The Carnot efficiency of the MHE is 


(28) 


Q. 


Tankl 


1034.21 

689.48 

344.74 

0 

-344.74 
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Fig. 8. Pressure variations in one pair of MHE heat exchangers (gauge pressure). 


Vldeal = 1 “ (29) 

where the values of Tc s and T Hs are the temperatures of the heat 
sink and heat source, respectively. 

The second law provides an indication of the irreversibilities 
within the system. The cycle exergy efficiency of the MHE based on 
the second law is defined as 




nth 

Vldeal 


(30) 


This can be calculated by substituting the cycle efficiency of the 
engine and Carnot efficiency into Eq. (30). The exergy efficiency 
provides an alternative interpretation of performance. When 
evaluating the energy efficiency, the same weight is assigned to 
energy whether it is electric, shaft work or low-temperature 
thermal energy. The exergy efficiency weighs the energy flows with 
respect to exergy [25]. Using the exergy efficiency, it can be better 
understood how internal and external irreversibilities affect the 
performance of the engine. The exergy efficiency (jje*; same as 
second law efficiency) for the MHE is defined as 


Vex = 


W e 

EXin 


(31) 


Exi n is the exergy entering the system. To calculate £x in , the 
exergy of hot water can be obtained from Eq. (12) (point 1 in Fig. 5). 
Thus, the value of Ex,„ is the exergy of the hot water stream that 
flows through the hot heat exchanger. Alternatively, Eq. (16) can be 
used to find £x !n . Considering Eq. (16), the denominator is smaller 
than the denominator of Eq. (28). Hence, the exergy efficiency of 
the system is higher than the system energy efficiency. This can be 
understood by low grade energy (such as waste heat) in the MHE 
that is converted to higher grade energy in the form of electricity. 


5. Results and discussion 

In this section, experimental data is obtained for the Marnoch 
heat engine, and compared against predictions from the thermo¬ 
dynamic model in the previous section. Fig. 8 shows how the 
pressure measurements vary within one pair of heat exchangers 
(hot and cold) during the operating period. This figure shows that 
the engine starts operation when the air pressure inside the first 
shell was 882.53 kPa, and the pressure of the other shell was 
103.42 kPa. Before starting the experiment, there was a 160 s delay. 
Within this period, as shown in Fig. 8, the pressure inside 


exchanger 1 dropped. The temperature of compressed air was 
initially higher than the ambient temperature, and the shell side of 
the heat exchanger was not perfectly insulated. Heat transfer from 
the air to the shell side, and from the shell side to the surroundings, 
occurred. In exchanger 2, the initial pressure inside the shell was 
close to the ambient pressure, so within the first 160 s, no 
considerable pressure change was observed in exchanger 2. The 
results in Fig. 8 refer to one operational cycle of the MHE and the 
measurement uncertainty (U) in the reported data is ±10 kPa. 

When the engine operation starts, due to heat transfer and air 
flow from exchanger 1 to exchanger 2, the pressure of exchanger 1 
decreases, and the pressure of exchanger 2 increases, simulta¬ 
neously. The curved sections in Fig. 8 show this period. When the 
piston cycle ends, the pressures inside the shells did not equalize. 
The remaining pressure difference was about 20.68 kPa in the 
experiment. This occurs because of friction losses in the piston 
assembly, so the remaining pressure difference could not move the 
piston. The pressure differences were determined from the exper¬ 
iments and predictions, to find the power output from the piston 
assembly and generator. The ideal gas law was used to predict the 
pressure differentials in the heat engine at various operating 
temperatures. 

Before each test, shells were filled with compressed air at 
different pressures. During the experiments, hot and cold water 
passed through the heat exchangers, in order to heat and cool them. 
Fig. 9 shows the total measured power output from the electric 
generator during one operation cycle. To obtain values for each 
pressure differential, the engine was operated at least three times 
for each initial pressure differential. The results shown in Fig. 9 are 
the average readings of one set of experiments. When the engine 
works under high pressure differences, the trend line in Fig. 9 
shows the electricity generated from the system increases. 

Fig. 10 shows the total mechanical power output from the piston 
assembly at specified pressure differences. The mechanical power 
outputs are determined by measuring the number of piston revo¬ 
lutions, within a certain time period, and the piston specifications 
from the manufacturer. The number of piston revolutions was 
found by measuring the number of changes (open/close) of the 
pneumatic valve. Comparing Figs. 9 and 10, the amount of power 
loss within the transmission system of the MHE can be observed. 

Fig. 11 shows the measured power output from the generator 
during one cycle of operation. The power is measured after each 
piston revolution [26]. The measured values for each curve in the 
graph represent one set of operating conditions of the heat engine. 
The parameter that has an effect on the measurement uncertainty 

































P. Saneipoor et al. / Energy 36 (2011) 6879-6889 


6887 



Fig. 9. Total power output from the generator at varying pressure differences. 
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Fig. 11. Measured power output with one exchanger pair at varying pressure 
differentials. 



is the instrumentation error. The measurement error of the 
instrument in the current setup is ±0.95 (W). The minimum 
pressure differential for these tests was 344.74 kPa, and the 
maximum was 1034.21 kPa. The power output from the generator 
shows fluctuations during the operation cycle, due to losses 
occurring in the flywheel. During operation at certain times, the 
speed of the flywheel is higher than the piston speed. Conse¬ 
quently, the energy output of the piston is not added to the 
flywheel energy. When the flywheel speed decreases, the power of 
the piston is added to the system again. 

In Fig. 12, the points represent the measured power output 
from the generator at two different pressure differentials. The 
continuous line shows the predicted power output from the 
generator at a 779.11 kPa initial pressure differential. The experi¬ 
mental pressure readings for this power prediction are presented 
in Fig. 8. To validate the predicted power output from the 


generator, the power output measurements from the upper and 
lower pressure differentials were used. This line is located 
between the power output readings of the initial pressure differ¬ 
ence of 689.48 and 1034.48 kPa. 

The predicted power output starts from a maximum value and 
drops over time gradually. However, the power measurements start 
from almost zero, increases higher, drops drastically and then the 
power reaches its maximum. The reason for this trend is the initial 
energy required to rotate the flywheel. This trend is visible in most 
of the MHE measured data (see Fig. 12). 

Fig. 13 shows a comparison of measured and predicted peak 
power outputs from the generator for varying pressure differen¬ 
tials. Seven observed peak power outputs at certain initial pres¬ 
sures are shown in this graph. The continuous solid line shows the 
predicted peak power outputs from the generator. The closest 



Fig. 10. Total power output from the piston assembly at varying pressure differentials. 



Fig. 12. Comparison of measured and predicted power outputs from the generator for 
varying pressure differentials during one operation cycle. 
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Fig. 13. Comparison of measured and predicted peak power outputs from the gener¬ 
ator for varying pressure differentials. 

match between the model and measured data occurs between the 
initial pressures of 600 and 700 kPa. At pressures around 1000 kPa, 
the model predicts higher power outputs from the system than the 
measured data. This occurs because the transmission efficiency was 
measured at 689.48 kPa. At higher pressure differentials, the rotary 
speed of the flywheel increases. Thus, more energy from the piston 
revolutions is not captured. At high pressure differentials, the 
transmission efficiency is less than 40%. 

Fig. 14 shows the Carnot efficiency of the system. The temper¬ 
ature ranges for the heat source and heat sink in this study are 
between 272 K and 373 K. In this range, water can add or remove 
heat from the system. In addition, in many areas of industry, waste 
heat in this temperature range is freely available. In Fig. 14, each line 



demonstrates the change of Carnot efficiency by changing the heat 
source temperature at a constant heat sink temperature. The main 
trend of the curves shows an increase in Carnot efficiency by 
elevating the source temperature. For instance, at a temperature of 
293 K and 298 K, for heat sink and source temperatures, the Carnot 
efficiency of a heat engine is 1.7%. When the temperature of the 
heat source is 353 K, the efficiency reaches 22.5%. The figure 
demonstrates that having the same temperature differences in 
lower temperature ranges gives a higher Carnot efficiency of the 
heat engine. When the heat sink temperature is 273 K, the effi¬ 
ciency of the heat engine is higher than 283 I< and 293 K. This trend 
is valid for temperature differentials below zero as well. Also, in 
Fig. 14, the maximum possible efficiency achievable from this 
system in temperature ranges of 273 K—373 K is 27.8%. Conse¬ 
quently, when water is used to add and remove heat from the 
system, the thermal efficiency of the system cannot exceed 27.8%. 

In the Marnoch heat engine, exergy losses occur in the system 
due to heat losses and mechanical component deficiencies. The 
electric power from the MHE is the exergy recovered by this heat 
engine. About 60% of the mechanical power output from the piston 
is lost due to friction in chains, belt, ball bearings and flywheel [26]. 
The existing mechanical loss in the MHE reduces the exergy output 
(recovered). In the current experimental apparatus, the heat loss 
from the heat exchangers is a significant portion of the total heat 
input to the heat exchangers [27]. When the ambient temperature 
is 298 K and the heat is supplied to the engine at a temperature of 
360 K, the dimensionless exergetic temperature factor is 0.17. It is 
assumed that the heat engine operates at an atmospheric pressure 
of 101.3 kPa. This indicates the exergy flow rate associated with the 
heat loss, £x/ oss , at atmospheric conditions is about 183.3 W. 

6. Conclusions 

This paper has presented, analyzed and performed experi¬ 
mental studies to demonstrate the promising performance of a new 
type of low-temperature heat engine, called a Marnoch heat 
engine. The heat engine is an energy conversion device that has the 
capability to recover low grade heat to electricity, using smaller 
temperature differentials than other existing methods. The MHE 
has strong potential to be used in many applications, such as waste 
heat recovery from thermal power plants, industrial facilities, 
furnaces and others. This paper presented new data for the thermal 
efficiency and power output of the MHE, as well as a predictive 
model that was validated against the experimental measurements. 
The studies indicate that higher pressure units will result in larger 
power output, higher efficiency and better economic viability of the 
system. 
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